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In this talk, we report an analysis of elliptic flow (^2), non-flow (^2), and V2 fluctuations 
{(Tv^) for the STAR collaboration at middle rapidity in y/s^ = 200 GeV Au-|-Au collisions. 
We compare the results to models of the fluctuations in the initial eccentricity of the 
overlap zone. 

Introduction: Understanding the early-time dynamics in heavy-ion collisions is an im- 
portant step in understanding the properties of the matter created in the collision overlap 
zone. Elliptic flow (^2) measurements are sensitive to the shape of the initial overlap zone [ 
[1] so f 2 fluctuations can reveal information about fluctuations and correlations in the ini- 
tial geometry. Those fluctuations and correlations will lead to a better understanding of 
the initial conditions of the collision evolution. We report on a study of V2 fluctuations 
based on 15 million events measured with the STAR detector [[2]- A similar analysis has 
been carried out by PHOBOS [|3] and results have been discussed in the literature [H]. 
Please note, data shown in this talk, showing small f 2 for events containing jets, was in 
error. 

Analysis method: The q-distribution method can be used to gain access to information 
about azimuthal correlations and f 2 fluctuations. The n*'* harmonic reduced flow vector 
is deflned as qn,x = -^^i cos^ncpi) and qn,y = -^Sj sin(n0j), where M is the number of 
tracks and 0, is the azimuth angle of a track with respect to a preferred emission axis [ 
[5]. The sums over the x and y components of the particle momentum are equivalent to 
random walks so that the central limit theorem is applicable for large M, in which case the 

distribution will be a 2-D Gaussian with widths that depend on particle correlations 
and V2 fluctuations. According to the definition of V2, if the reaction-plane is the preferred 



axis, the 2-D Gaussian will be shifted along x by ^/Mvn with widths: ~ | (1 + M5n) 
and ay ~ | (1 + 5„ = (cos(nA0)) — are correlations unrelated to the reaction 

plane [ [5j . Neither the preferred axis, nor the reaction plane direction are known on an 
event-by-event basis, so we calculate the magnitude The 2-D Gaussian then becomes: 



dN 
qndqv 




(1) 



where P and Ik are gamma and modified Bessel's functions. If we assume the preferred 
axis for the dN/d\q\ is the reaction-plane, then V2{q} = V2 and V2 fiuctuations can be 
introduced by convoluting the above distribution with a V2 distribution {dN/dv2)- This 
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introduces an added width to dn/d\q2\: ~ i + + '^Ma'^^j, where a^^ is the 
r.m.s. width of dN/dv2- 62 and a^^ contribute to the width with a factor of M, so it is not 
possible to extract these quantities independently as previously thought [i6j|. Rather, we 
define the dynamic width of dN/d\q2\ as cr^y„ = 62 + Scx^^' which can be readily extracted 
from data. 

We don't know a priori what the preferred emission axis is. For peripheral Au+Au or 
p+p collisions, the preferred axis may be predominantly determined by the thrust axis 
of outgoing partons. However, for mid-central and central collisions, the four-particle 
cumulant f2{4} (a higher multiplicity correlation analysis [[?]) decreases with eccentricity 
suggesting the geometry of the overlap determines the preferred axis. 

Distinguishing between cr^^ and 62 requires knowledge of the reaction-plane or infor- 
mation about the A0-, At]-, charge-sign- or multiplicity-dependence of non-flow. STAR 
addresses this in two general ways, 1) analyzing higher multiplicity correlations as in a 
q-distribution or higher order cumulant analysis [ O [1] or 2) analyzing the A0, At], or 
charge-sign dependence of correlations. Examples of 2) include measuring two-particle 
correlations across a wide r^-gap [ [8] , correlating only like-sign particles, or by fitting the 
A0, Ar] correlation density to isolate the cos2(A0) term V2{2D} [[9]. But the small A(j) 
correlations in Au+Au collisions stretch over a large Arj range [[10] and it's possible that 
the space-momentum correlation inducing mechanism responsible for f 2 may contribute 
to the near-side correlation by translating quantum- or geometry-fluctuations from the 
initial overlap to correlations in momentum-space [ |TTl |T2] linking the ridge and V2 [ US] • 

For dN/dv2, we consider a Gaussian and Bessel-Gaussian {BG{v bg, '^vbg)) [tlSl- Convo- 
luted with Eq. [H both give an equally good description of dN/d\q\. The mean and r.m.s. 
of the EG or the parameters vbg and 

c^vgQ can be reported. We find that the BG mean 
and r.m.s. are within errors the same as the Gaussian mean and r.m.s. The use of a BG 
shape for the dN/ dv2 is motivated by the shape of the participant eccentricity [ dis- 
tribution in a Glauber Monte Carlo model In the case that V2 (x. 62, the parameters 
extracted from the BG are either related to the participant-plane or the reaction-plane. 
The parameter vbg represents the V2 computed with respect to the reaction-plane, while 
{V2) represents V2 computed with respect to the participant-plane [[16]. Here, we report 
the mean and rms width of a Gaussian and compare it to models of Epart- 

Results: Fig. [T] (left) shows our fits to dN/d\q2\ where \q2\ is calculated from three track 
samples: full-acceptance (qjfull}), reduced-acceptance {q{ri/2}), and full-acceptance like- 
sign (q{ like-sign}). The different samples are expected to contain different contributions 
of correlations; HBT, resonance decays, and parton fragments, which contribute to non- 
zero values of 62- The inset shows the ratio of each distribution to (iA^/(i|q{full}|. The 
widths are ordered as follows: adyn{<i{v/'^}} > crdj/n{ci{full}} > o"rfy„{g{like — sign}}. This 
is expected because the like-sign sample removes most correlations caused by resonance 
decay. Jet fragments exhibit charge-ordering so those correlations are also reduced. The 
ri/2 sample gives a greater weight to the short-range correlations so 62 is enhanced. 

Fig. [1] (right) shows the dynamic width of the full-acceptance dN/d\q\. The tan band 
includes systematic uncertainties from acceptance, finite multiplicity and comparisons of 
methods including the dynamic width from the cumulants f2{2}^ — f2{4}^ = ^2 + 2o"^2- 
It was shown in ref. [[T6] that within the precsion of our measurements we cannot extract 
more information from dN/d\q\ than that already found from the cumulants. In particular, 
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Figure 1. Left: dN/d\q\ for like-sign tracks (q{ like-sign}), tracks within one hemisphere 
of the STAR acceptance {q{r]/2}), and all tracks with 0.15 < pt < 2.0 GeV/c (q{full}). 
The inset shows the ratio of each distribution to q{full}. Right: The dynamic width 
of dN/d\q\ with width from cumulant analysis. f2{2}^ — V2{2D}'^ gives the non-flow S2 
inferred from a fit to the 2D correlation densities [i9j. 

if the f 2 fluctuations are Gaussian, then the equality V2{2}^ — V2{4}^ = 62 + 2al^ is exact. 
Published STAR f2{2} and f2{4} values vs. centrality have been calculated by weighting 
each event by the number of tracks in the event [ [1]. Those are not comparable to 
the q-distribution analysis where each event is weighted equally. We therefore compare to 
cumulant data without the multiplicity weighting and find good agreement. The minimum 
value of S2 shown in Fig. [T] (right) is calculated from the difference between {q*q) calculated 
with same-charge pairs and unlike-sign pairs. We also show 62{2D} = W2{2}^ — V2{2D}^ 
which shows that a large component of a'^y^ is accounted for by the near-side peak in the 
2D correlations [[9]. Fig. [2] (left) shows the allowed values for (t>2) and ay^. The case of 
zero fluctuations is not excluded and an upper limit is reported. 




impact parameter <b> (fm) impact parameter (b> (fm) 

Figure 2. Left: The tan band shows the allowed values of (^2) from this analysis and the 
upper limit on a^^- Right: Upper limit on (Jy^/{v2) compared to models of eccentricity. 

Fig. [2] (right) shows a^.^/ {^2) calculated using a Gaussian for dN/dv2. The upper limit on 
this ratio is compared to several Monte Carlo models of Epart fluctuations. In the nucleon 
Glauber models, the eccentricity e = {y"^ — x"^) / {y"^ + x^) is calculated based on the x 
and y cooridinates of participant nucleons. For the constituent quark Glauber model, 
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each nucleoli is treated as three constituent quarks bound within the nucleon, and the 
eccentricity is calculated based on the constituent quarks that participate in the collision. 
All models lie within the allowed range while the CGC model [|T7] and the model based 
on constituent quarks both have smaller relative widths. The nucleon participant model 
leaves little room for other sources of fluctuations and correlations beyond the initial 
geometric ones. The large near-side peak observed in two-particle correlations contradict 
the idea that all or most of the width of a'^y^ is dominated by a^^ suggesting that the CGC 
or constituent-quark model may be preferred. Recently it was proposed that correlations 
and fluctuations in the initial conditions may also contribute to the near-side ridge [[T2]. 

Conclusions: The dynamic width of dN/dq has been presented. The width reflects 
correlations and f 2 fluctuations. An upper limit on a^^ has been presented and compared 
to models of Spart fluctuations. From our comparison to the Glauber MC we conclude 
that either the density fluctuations in the initial state are manifested in both the near- 
side peak and in a^,^, or the initial overlap density is smoother than would be expected 
from a Glauber model. 
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